Data-driven identification of predictive risk biomarkers for subgroups of osteoarthritis using interpretable machine learning.

Osteoarthritis (OA) is increasing in prevalence and has a severe impact on patients' lives. However, our understanding of biomarkers driving OA risk remains limited. We developed a model predicting the five-year risk of OA diagnosis, integrating retrospective clinical, lifestyle and biomarker data from the UK Biobank (19,120 patients with OA, ROC-AUC: 0.72, 95%CI (0.71-0.73)). Higher age, BMI and prescription of non-steroidal anti-inflammatory drugs contributed most to increased OA risk prediction ahead of diagnosis. We identified 14 subgroups of OA risk profiles. These subgroups were validated in an independent set of patients evaluating the 11-year OA risk, with 88% of patients being uniquely assigned to one of the 14 subgroups. Individual OA risk profiles were characterised by personalised biomarkers. Omics integration demonstrated the predictive importance of key OA genes and pathways (e.g., GDF5 and TGF-ß signalling) and OA-specific biomarkers (e.g., CRTAC1 and COL9A1). In summary, this work identifies opportunities for personalised OA prevention and insights into its underlying pathogenesis.

A polygenic biomarker to identify patients with severe hypercholesterolemia of polygenic origin.

BACKGROUND: Severe hypercholesterolemia (HC, LDL-C > 4.9 mmol/L) affects over 30 million people worldwide. In this study, we validated a new polygenic risk score (PRS) for LDL-C. METHODS: Summary statistics from the Global Lipid Genome Consortium and genotype data from two large populations were used. RESULTS: A 36-SNP PRS was generated using data for 2,197 white Americans. In a replication cohort of 4,787 Finns, the PRS was strongly associated with the LDL-C trait and explained 8% of its variability (p = 10-41 ). After risk categorization, the risk of having HC was higher in the high- versus low-risk group (RR = 4.17, p < 1 × 10-7 ). Compared to a 12-SNP LDL-C raising score (currently used in the United Kingdom), the PRS explained more LDL-C variability (8% vs. 6%). Among Finns with severe HC, 53% (66/124) versus 44% (55/124) were classified as high risk by the PRS and LDL-C raising score, respectively. Moreover, 54% of individuals with severe HC defined as low risk by the LDL-C raising score were reclassified to intermediate or high risk by the new PRS. CONCLUSION: The new PRS has a better predictive role in identifying HC of polygenic origin compared to the currently available method and can better stratify patients into diagnostic and therapeutic algorithms.

Identification of disease-associated loci using machine learning for genotype and network data integration.

MOTIVATION: Integration of different omics data could markedly help to identify biological signatures, understand the missing heritability of complex diseases and ultimately achieve personalized medicine. Standard regression models used in Genome-Wide Association Studies (GWAS) identify loci with a strong effect size, whereas GWAS meta-analyses are often needed to capture weak loci contributing to the missing heritability. Development of novel machine learning algorithms for merging genotype data with other omics data is highly needed as it could enhance the prioritization of weak loci. RESULTS: We developed cNMTF (corrected non-negative matrix tri-factorization), an integrative algorithm based on clustering techniques of biological data. This method assesses the inter-relatedness between genotypes, phenotypes, the damaging effect of the variants and gene networks in order to identify loci-trait associations. cNMTF was used to prioritize genes associated with lipid traits in two population cohorts. We replicated 129 genes reported in GWAS world-wide and provided evidence that supports 85% of our findings (226 out of 265 genes), including recent associations in literature (NLGN1), regulators of lipid metabolism (DAB1) and pleiotropic genes for lipid traits (CARM1). Moreover, cNMTF performed efficiently against strong population structures by accounting for the individuals' ancestry. As the method is flexible in the incorporation of diverse omics data sources, it can be easily adapted to the user's research needs. AVAILABILITY AND IMPLEMENTATION: An R package (cnmtf) is available at https://lgl15.github.io/cnmtf_web/index.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Properties of human genes guided by their enrichment in rare and common variants.

We analyzed 563,099 common (minor allele frequency, MAF≥0.01) and rare (MAF < 0.01) genetic variants annotated in ExAC and UniProt and 26,884 disease-causing variants from ClinVar and UniProt occurring in the coding region of 17,975 human protein-coding genes. Three novel sets of genes were identified: those enriched in rare variants (n = 32 genes), in common variants (n = 282 genes), and in disease-causing variants (n = 800 genes). Genes enriched in rare variants have far greater similarities in terms of biological and network properties to genes enriched in disease-causing variants, than to genes enriched in common variants. However, in half of the genes enriched in rare variants (AOC2, MAMDC4, ANKHD1, CDC42BPB, SPAG5, TRRAP, TANC2, IQCH, USP54, SRRM2, DOPEY2, and PITPNM1), no disease-causing variants have been identified in major, publicly available databases. Thus, genetic variants in these genes are strong candidates for disease and their identification, as part of sequencing studies, should prompt further in vitro analyses.

Serum Adipsin Levels throughout Normal Pregnancy and Preeclampsia.

Adipsin is a protease produced at high levels by adipose tissue. It is involved in complement activation and metabolic control. The objective of this study was to determine the changes in adipsin levels during different stages of normal pregnancy, and its association with obstetric outcomes, such as preeclampsia. This nested case-control study in a longitudinal cohort included normal pregnant (n = 54) and preeclamptic (n = 18) women, both followed throughout pregnancy. Additionally, some of the normal pregnant women were followed up three months postpartum (n = 18). Healthy non-pregnant women were also studied during their menstrual cycle (n = 20). The results of this study show that in healthy non-pregnant women, adipsin levels did not change significantly during the menstrual cycle. In normal pregnant women, adipsin levels were lower (p < 0.01) when compared with non-pregnant healthy women, but these serum levels increased again during postpartum (p < 0.001). Adipsin levels were significantly elevated in preeclamptic women in late pregnancy (P < 0.01). A significant correlation was not found between leptin and adipsin during the three periods of gestation studied in healthy pregnant and preeclamptic women. Our results suggest that adipsin may be involved in pregnancy-associated metabolic changes. Moreover, the increase of adipsin levels towards late gestation in preeclamptic women could be related to the pathophysiology of this disease.

Maternal serum omentin-1 profile is similar in humans and in the rat animal model.

Omentin-1 is an adipocytokine with anti-inflammatory activity that has been associated with different metabolic disorders. The aim of this study is to investigate the serum profiles of omentin-1 throughout human and rat pregnancy. Serum omentin-1 levels were determined by ELISA in a prospective cohort study of healthy pregnant women (n=40) during the three trimesters of pregnancy and in twenty healthy non-pregnant women during the follicular and luteal phase of the menstrual cycle. In addition, serum omentin-1 levels were measured in rats during different periods of pregnancy (gestational days 8, 12, 16, 19, and 21) and in an age-matched control (virgin) group of rats (n=12rats/group). Finally, immunohistochemistry was used to demonstrate the presence of omentin-1 protein in human and rat placenta. Omentin-1 immunoreactivity was detected in cytotrophoblasts, syncytiotrophoblasts, sparse Hofbauer cells, and endothelial cells of the stem villi of human placenta. Additionally, it was detected in the labyrinthine trophoblast and yolk sac layer of the rat placenta. Human and rat serum omentin-1 levels were significantly lower in the late gestational period when compared with the non-pregnant women and virgin rats (p<0.05). Serum omentin-1 changes were not significant throughout the gestation in both species (p>0.05). Human serum omentin-1 levels have an inverse relationship with triglyceride levels during pregnancy. Our findings have not determined the exact role of omentin-1 during pregnancy, concerning the metabolic control of triglycerides and other energy sources. Whether omentin-1 decrease implies a regulatory function is still not clear. Further studies are needed to address this issue and determine the role of omentin-1 in metabolic adaptations during normal human and rat pregnancy.

Maternal Serum Meteorin Levels and the Risk of Preeclampsia.

BACKGROUND: Meteorin (METRN) is a recently described neutrophic factor with angiogenic properties. This is a nested case-control study in a longitudinal cohort study that describes the serum profile of METRN during different periods of gestation in healthy and preeclamptic pregnant women. Moreover, we explore the possible application of METRN as a biomarker. METHODS AND FINDINGS: Serum METRN was measured by ELISA in a longitudinal prospective cohort study in 37 healthy pregnant women, 16 mild preeclamptic women, and 20 healthy non-pregnant women during the menstrual cycle with the aim of assessing serum METRN levels and its correlations with other metabolic parameters. Immunostaining for METRN protein was performed in placenta. A multivariate logistic regression model was proposed and a classifier model was formulated for predicting preeclampsia in early and middle pregnancy. The performance in classification was evaluated using measures such as sensitivity, specificity, and the receiver operating characteristic (ROC) curve. In healthy pregnant women, serum METRN levels were significantly elevated in early pregnancy compared to middle and late pregnancy. METRN levels are significantly lower only in early pregnancy in preeclamptic women when compared to healthy pregnant women. Decision trees that did not include METRN levels in the first trimester had a reduced sensitivity of 56% in the detection of preeclamptic women, compared to a sensitivity of 69% when METRN was included. CONCLUSIONS: The joint measurements of circulating METRN levels in the first trimester and systolic blood pressure and weight in the second trimester significantly increase the probabilities of predicting preeclampsia.

Longitudinal analysis of maternal serum Follistatin concentration in normal pregnancy and preeclampsia.

OBJECTIVE: Follistatin (FST) is a regulator of the biological activity of activin A (Act A), binding and blocking it, which could contribute to the modulation of its pro-inflammatory activity during pregnancy. We sought to investigate, in this nested case-control study, FST serum levels during normal pregnancy and correlate it with the FST profile in preeclamptic pregnant women, normal pregnant women followed 3 months postpartum and eumenorrheic nonpregnant women throughout the menstrual cycle. SUBJECTS AND METHODS: Follistatin serum levels determined by ELISA, biochemical and anthropometric variables were measured in normal pregnant (n = 28) and preeclamptic (n = 20) women during three periods of gestation. In addition, FST serum levels were measured in a subset of normal pregnant women (n = 13) followed 3 months postpartum and in eumenorrheic nonpregnant women (n = 20) during the follicular and luteal phases of the menstrual cycle. RESULTS: Follistatin serum levels in the eumenorrheic nonpregnant and postpartum group were significantly lower when compared to levels throughout gestation (P < 0·01). Serum FST levels increased in each period of pregnancy analysed, being significantly higher towards the end of gestation (P < 0·01). FST levels were lower in late pregnancy in preeclamptic women compared to normal pregnant women (P < 0·05). Finally, FST levels were higher in the luteal phase when compared with the follicular phase of the menstrual cycle (P < 0·05). CONCLUSIONS: These analyses would permit the consideration that changes in FST levels during pregnancy contribute to the control of the Act A system.